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ABSTRACT 

The  Mia  results  of  a  structure!  Analysis  program  concerning 
a  conventional  submarine  pressure  hull  configuration  consisting  of  a 
rib-stiffened  cylindrical  shell  with  a  reinforced  circular  penetration 
under  hydrostatic  pressure  are  discussed.  The  mathematical  solution 
in  series  form  has  been  obtained  by  superposing  the  analytical  solutions 
for  (a)  a  long,  unatlffeoed,  end  unperforated  circular  cylindrical 
thin  shell,  closed  at  the  end*,  under  external  hydrostatic  pressure, 

(b)  a  long,  unatiffened,  and  imperforated  circular  cylindrical  thin  shell 
under  a  prescribed  number  of  arbitrary  radial  line  loads,  and  (c)  a 
long,  unstiffened,  circular  cylindrical  shallow  thin  shell  under 
arbitrary  loading  along  the  boundary  of  a  circular  penetration.  During 
the  computation,  following  truncation  of  the  series,  the  method  of 
least  squares  la  employed  in  solving  for  the  integration  constants 
determined  by  the  boundary  conditions  preocrlbed  along  the  ribs  and 
the  reinforced  penetration.  The  analysis  has  been  coded  and  the 
numerical  results  generated  by  the  use  of  a  digital  computer  for  an 
unatiffened  as  well  as  a  rib-stiffened  shell  with  a  reinforced  penetration 
are  presented  and  compared  graphically  with  experimental  data  available 
from  certain  photoelastic  model  tests.  The  corresponding  zone  of 
influence  of  the  penetration  and  the  state  of  stress  concentration 
about  it  are  delineated. 
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*  subscript  denoting  solution  for  a  plain,  closed  cylindrical 
shall  under  pressure. 

*  coefficients  depending  on  v,  ic,  m;  i  *  1,  2,  ...  ,  5 

■  roots  of  A  in  lover  half  of  tooplex  plane;  s  ■  1,  2,  3,  4 

(o) 

*  complex  conjugates  of 

»  area  of  rib  cross  section 
»  area  of  ring  cross  section 

■  integration  constants 

*  integration  constants 

»  subscript  denoting  solution  for  &  plain  cylindrical  shall 
under  multiple  radial  line  loads 

*  coefficients  depending  on  v>,  xt  m;  1*1,  2,  ...  ,  9 
c  constant  defined  in  Kq.(2  6> 

*  quantities  defined  in  Eq.(35) 

«  inverse  of  B  .. 

B'jjL 

»  subscript  denoting  solution  for  a  plain  cylindrical  shell 
loaded  along  a  circular  penetraticu 

*  rigid  body  displacement  of  i-th  pair  of  ribs 

»  torsional  rigidity  of  ring 

*»  functions  defined  in  Eq.(3j) 

»  functions  defined  in  Eq. (3l) 
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•  area  moment  of  Inertia  of  rib  croft*  section  relative  to 
centroidal  principal  axis  normal  to  plane  of  rib 

*  area  aoucut  of  inertia  of  ring  cross  section  relative  to 

1- 1  axis,  Fig. 3 

•  area  moment  of  inertia  of  ring  cross  section  relative  to 

2- 2  axis,  Fig. 3 

•  imaginary  part 

*  subscript  integer 
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»  subscript  integer 


*  half  distance  between  j-th  pair  of  ribs;  j  -  1,  2,  ,  N 

j  -  1,  Fig. I 


“  subscript  integer 
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M0n  **  shei)L  couple  resultants  in  r»,  8  coordinate  system.  Fig. 2 
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shall  couple  resultants  in  p,  coordinate  system.  Fig. 2 
subscript  integer 

number  of  pairs  of  ribs  located  symmetrically  across  trana 
verse  axle  of  shell 
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i.  iraroucnoif 

Current  interest  in  the  structural  analysis  of  deep  diving  Manned 
submsrslbles  end  of  submerged  habitats,  such  as  naval  or  commsrdal 
submarines  end  installations  and  oceanographic  research  vehicles  and 
laboratories,  is  due  to  evolving  tactical  concepts  and  easrging  eco- 
nosdLc  and  scientific  goals  in  the  development  of  underwater  resources 
and  the  exploration  of  the  underveter  environment.  It  is  generally 
recognised  thet  the  presence  of  structural  penetrations,  such  as 
missile  tubes,  hatches,  locks,  observation  ports,  valves,  and  various 
other  functional  connections,  constitutes  a  source  of  weakness  in  any 
pressure  hull.  Related  problems  arise  in  the  pressure  vessel  and 
pipeline  technologies  where  structural  Intersections  are  common.  A 
comprehensive  review*  of  past  research  on  the  subject  Is  available. 

In  recent  years,  considerable  analytical  and  experimental  data 

on  the  stress  analysis  of  pressure  vessels  with  intersections  have 
2—4 

been  obtained  in  support  of  programs  for  the  development  of  asso¬ 
ciated  design  codes.  The  results  make  possible  the  planning  of  further 
experiments  and  the  partial  verification  of  information  obtained  by 
sore  powerful  methods^  of  computational  stress  analysis. 

Pressure  hulls  with  penetrations  have  also  been  Investigated 
by  both  analytical  and  experimental  means  in  the  course  of  recent 
programs  of  submarine  structural  research  involving  conventional 
hull  configurations.  The  present  vosdc,  part  of  such  efforts,  concerns 
the  analysis  of  a  rib-stiffened  cylindrical  shell  with  a  reinforced 
circular  penetration  subjected  to  hydrostatic  pressure.  The  special 
caass  of  a  rib-stiffened  shell  with  an  unreinforced  penetration^. 
Including  the  presence  of  multiple  stiffeners7,  and  that  of  an  unstlf- 

g 

fened  shell  with  a  reinforced  penetration  have  been  reported  earlier. 

q 

Other  numerical  results  on  the  last  problem  are  also  available.  Experi¬ 
mentally,  all  of  the  above  pressure  hull  configurations  have  been 
•tudled.10-12 
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The  configuration  considered  here,  as  shown  In  Fig.l,  cons  lata 
of  a  long,  circular  cylindrical,  thin  shall  of  uniform  trail  thickness 
t  and  Biddle  surface  radius  R,  dosed  at  the  ends,  with  several 
identical,  internal,  circumferential  stiffeners  of  centroidal  radius 
R9 ,  and  a  single,  centrally  located,  reinforced  circular  penetration 
of  outer  radius  po,  closed  by  a  flat  nenbrane.  The  loading  consists 
of  external  hydrostatic  pressure  p,  and  the  deformation  is  taken  snail 
and  static.  The  dineneionlese  axial  and  circumferential  coordinates, 
n  end  0,  respectively,  are  centered  at  the  point  of  intersection  0 
of  the  longitudinal  end  transverse  axes  of  the  shell,  as  shown  in 
Fig.l,  with  R»j  a  measure  of  length.  The  corresponding  axial,  circum¬ 
ferential,  and  radial  displacements ,  U,  V,  and  W,  respectively ,  stress 
wulUnf  V  He,  Nne,  V  Q^,  Qg,  couple  resultants  *V  Me’  Mne* 

end  an  arbitrary  radial  load  p1  for  an  element  of  the  shell  middle 
surface  are  indicated  in  Fig.  2  where  the  intended  sense  of  a  couple 

is  also  given.  The  surface  coordinates  p  and  +,  representing  polar 

2 

coordinates  defined  in  the  plane  obtained  by  developing  the  cylindrical 
middle  surface,  ere  centered  at  the  point  0*  directly  above  point  0, 
at  the  intersection  of  the  main  generator  and  the  transverse  axle  of 
the  shell,  as  shown  In  Fig.l,  with 

Hi)  ■  pcoe  RO  *  pain  (1) 

The  corresponding  displacements  U  ,  U.,  stress  resultants  N  ,  N.,  N  , 

P  7  P  $  Pt 

N^,  Qp,  Q^,  and  couple  resultants  H^,  M^,  for  an  element  of 

the  shell  middle  surface  are  indicated  in  Fig. 2.  Partial  derivatives 
will  be  denoted  by  a  comma  followed  by  the  particular  coordinates  as 
subscripts. 

The  circumferential  stiffeners  or  ribs  are  assumed  to  act  as 
slender,  curved  beams  possessing  a  plane  of  longitudinal  symmetry  and 
undergoing  only  flexural  and  extensions!  deformations  and  radial  inter¬ 
action  with  the  shell.  The  relevant  geometrical  properties  of  a  normal 
cross  section  of  a  rib  are  then  the  area  A  and  the  area  moment  of 
inertia  I  with  respect  to  the  centroidal  principal  axis  normal  to  the 
plane  of  the  rib. 
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A  representative  reinforced  circular  penetration  or  ring 
intersecting  the  shell  is  shown  in  Fig* 3  in  the  transverse  plane 
of  intersection  normal  to  the  aide  of  the  shell*  The  ring  is 
assumed  to  resist  flexural,  extensions!,  and  torsional  deformations* 

A  normal  cross  section  of  the  ring  has  an  area  A  with  an  axis  of 
symmetry  2-2  in  the  plane  of  the  ring.  The  axis  1-1  is  the  other 
centroldal  principal  axis  of  the  crocs  section,  located  a  distance 
p  from  the  transverse  axis  of  the  shell*  The  area  moments  of 
inertia  of  the  cross  section  with  respect  to  the  1-1  and  2-2 
axes  are  denoted  hy  ^  and  Ig,  respectively*  The  distance  e  from 
the  axis  1-1  to  the  point  of  intersection  of  the  ring  and  the 
shell  is  assumed  to  be  independent  of  t*  The  distance  e  is 
measured  from  the  axis  1-1  to  the  point  of  Intersection  of  the 
ring  and  the  flat  membrane  covering  the  opening.  The  vertical 
distance  from  the  plane  of  the  ring  to  points  above  It  on  the 
shell  middle  surface  at  p  ■  pQ  is  denoted  by  h(+),  and  the  depth 
of  the  ring,  by  h^« 

The  materials  of  the  shell,  ribs,  and  ring,  assumed  different 
in  the  analysis,  are  homogeneous,  isotropic,  and  elastic.  The 
material  constants  are  the  Young's  moduli  E,  E' ,  end  E,  and  Poisson's 
ratios  v,  v' ,  and  *  for  the  shell,  rib,  and  ring,  respectively* 

The  torsional  rigidity  of  the  ring  depends  on  E,  v,  end  the  geometry 
of  the  cross  section,  and  Is  denoted  by  C* 

The  solution  for  the  total  state  of  deformation  of  the  shell 
is  obtained  by  superposing  the  intermediate  solutions  for  (a)  s  long, 
unstlffened,  and  unparf orated  circular  cylindrical  thin  shell,  closed 
at  the  ends,  under  external  hydrostatic  pressure,  (b)  a  long,  unstif¬ 
fened,  and  unperforated  circular  cylindrical  thin  shell  under  a 
prescribed  number  of  arbitrary  radial  line  loads,  such  as  those 
exerted  by  the  ribs,  and  (c)  a  long,  unstlffened,  circular  cylindrical 
shallow  thin  shell  under  arbitrary  leading  along  the  boundary  of 
a  circular  penetration.  Quantities  referring  to  these  intermediate 
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solutions  (a),  (b),  and  (c)  arm  denoted  by  the  subscripts  a,  b9 
and  c,  respectively*  Shall  quantities  without  auch  subscripts 
refer  to  the  solution  for  the  total  state  of  deformation  evaluated, 
in  general,  by  summing  results  from  cases  (a),  (b),  and  (c) 
above.  This  total  solution  for  the  shall  swat  fulfill  specified 
boundary  conditions  along  the  riba  and  the  reinforced  penetration. 

Tor  instance,  the  loads  exerted  by  the  shell  on  the  ring  are 

«• 

shewn  in  Fig.  3  where  and  N^  are  respectively  the  effective 
transverse  and  membrane  shear  stress  resultants  defined,  in 
general,  as 


% 


s  %  ~ 


P_1  H 


PM 


N 


P* 


SHP*  +  R 


-X 


V°8 


with 


(2) 


K 


-X 


P>P 


<V 


w*> 


(3) 


The  final  form  of  the  solution  is  expressed  in  terms  of  Fourier 
series  expansions  in  the  6  or  +  coordinates,  with  coefficients  that 
are  functions  of  the  n  or  p  coordinates,  respectively,  end  contain 
arbitrary  integration  constants  and  various  constant  parameters. 

Daring  the  computation,  the  different  series  are  truncated  and  the 
method  of  least  squares  is  employed  in  solving  for  the  appropriate 
number  of  integration  constants  from  the  boundary  conditions.  Con¬ 
vergence  la  verified  numerically. 

The  numerical  data  reported  here  refer  to  the  state  of  stress 
in  an  unstlffened  and  a  rib-stiffened  cylindrical  shell,  respectively, 
in  the  presence  of  a  reinforced  circular  penetration.  The  results 
are  also  utilized  In  calculating  the  state  of  stress  concentration 
in  the  shell.  For  this  purpose,  a  dimensionless  principal  9tre&£ 
concentration  factor  N  is  defined,  in  general,  £? 


H  s  nax  |  ox,  a2  |  /ms  j  or  o2  j  , 


(4) 
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where  a^9  cr2  are  the  principal  stresses  at  a  point  of  tha  ahall 
with  tha  riba  and  tha  panatration  present  ,  while  3p  are  tha 
principal  atraaaaa  at  tha  asms  point  in  tha  absence  of  tha  riba 
and  tha  penetration*  Tha  former  are  given  by 

2*1,2  "  «p  +  %  4  “V  +  4  '*♦  ^  ’  (5> 

in  terms  of  the  normal  and  ahearing  stresses  ap,  o^,  and  *p^, 
respectively,  in  the  p,  4  coordinate  ays  ten,  while  the  latter,  for 
the  loading  considered  here,  correspond  to 


■“  I  «2  I  ■  pR/t  * 


(6) 


Denoting  by  a  the  thickness  coordinate  of  the  shell,  measured 
positive  outward  from  the  middle  surface,  with 


N  2  o  t  •»  N  -  12rM  /t  , 

P  P  P  P 

S  of  -  -  12*M^/t2, 

s» 1  V "  V  *  la,V‘!  • 


it  follows  dram  Eqo.(4)-(6)  that 

5  *  (2pR)~*  mas  \  Sp  4  ±  [  <Np  -  N^)2  +  4N2^  ]** 


(7) 


(8) 


From  Eq* (8) ,  setting  *  equal  to  *st,  -*gt,  aud  zero  in  Eq.(7),  the 
particular  values  of  5  at  the  outer,  inner,  and  middle  surfaces  of 
the  shell,  respectively,  are  determined . 


2.  CLOSED  CYLINDRICAL  SHELL  UNDER  HYDROSTATIC  PRESSURE 

For  a  long,  unstiffened,  and  imperforated  circular  cylindrical 
thin  shell,  closed  at  the  ends,  under  external  hydrostatic  pressure, 
the  solution  is 
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%•*,  \6  -  -pR,  o,  -  (V  -  J»)  pR2(Et)'ln  , 


Ha  -  *S(v  -  2)  pR2{Bt)"1 
In  the  tip  0  coordinate  system,  and 


M*p  *  (co*  2*  "  3)  *  (3  +  cos  2*)  , 

%♦  ■  "  -1»R  sln  *♦  . 
u«o  * (v  ■  ^  p4®6)"1**  c<*2  ♦ . 

U«*  “  ^  pR(Bt)"\>  2*  , 

Wa  -  5s(v.  -  2)  pR2  (Et)_1 

J 

in  the  p,  +  cor-dlnate  system,  all  other  quantities  being  equal 
to  zero,  so  that  from  Eqs. (2)  and  (3) , 


o  *  o,  n  *  n 

ap  •  ap(f  apf 


3.  CYLINDRICAL  SHELL  UNDER  MULTIPLE  ARBITRARY  RADIAL  LINE  LOADS 


Por  a  long,  unstiffened,  and  unperforated  circular  cylindrical 

thin  shell  under  an  arbitrary  radial  load,  the  governing  displacement 

13 

equations  of  equilibrium  are 


Vn*  +  **(1  ‘  v)  Vee  +  35(1  +  v)  vb.ev,  +  *b,n 

+<cft(l  -  >0  ob>96  -  Hb>ni)n  +  Wi  -  v)  wb>n9f }  -  o  , 

*a  +  v)  Vne  +  vb,ee  +  «l  "  v)  +  V 

+  4“  -  v)  vb.nn  '  *3  -  v)  Xme]  ”  0  * 
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v0b,n  +  Vb.e  +  «b  +  -  V)  0^,  - 

~  HO  ~  v)  t i^r)6  +  ^,6609  +  ®^>,nn99  +  ,nnwi 

+  ^.ea +  V  "  p'r2/I)  * 

k  5  im2  , 

where 


D  S  it/a  -  v2)  ,  K  =  Et3/12(l  -  V2) 


(13) 


are  respectively  the  extenslonal  and  flexural  rigidities  of  the  shell. 

In  the  particular  case  of  an  arbitrary  radial  line  load, 
solutions  of  Eq.(12)  decaying  to  zero  as  n  •>  *•  are  sought. 

Bxpmdiag  the  line  load  exerted  on  the  shell  by  a  single  rib 
locited  at  n  ■  0,  and  the  displacements  in  Eq.(12),  in  terms  of 
suitable  Fourier  series  in  6, 


<?'>  V  V  *  Lo  {Pm5(n>*  Vn)*  Vn)  '  008  °e . 
Vb  ‘  Lo  VaM  #in  “9  * 


(14) 


vaere  m  i  an  integer,  p^  are  arbitrary  constants  denoting  inter¬ 
action  iojd  amplitudes,  $(n)  la  the  Dirac  delta  function,  and  Uffl, 
V  ,  an  displacement  amplitudes,  and  applying  to  Eq.(12)  the 
Fourier  transform  defined,  for  an  arbitrary  function  f(n),  as 


f  (a)  « 


m 

t««i 


f(n)  exp  ian  dn 


(15) 


where  a  to  the  transform  variable  corresponding  to  n,  and 
i  s  (-1)  *.  a  system  of  equations  for  UQ(a)  s  Vm<a)t  and  Wm(a) 
is  obtain  d  the  solution  of  which  Is**  of  the  form 

i\,  Vo,  wffi>  «  pjMf1#  f2,  f3}/(2<rDA  ,  (16) 
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liMrt 

A  (a)  5  (a®  +  a^a*  +  +  tjO2  +  a^)  , 

jjtojs  iaO^a4  +  b2a2  +  b3)  , 

f2<«>  s  \a*  +  b5o2  +  b6  , 

f3(«)  s  b7«'  +  b8«2  +  b,  , 

and  the  coefficient a  a^,  1  «  1,  2,  ,  5,  and  b^,  1  *  1,  2,  ...  ,9, 

depend  on  vt  k,  b,  It  should  be  noted  that  p^  s  0  if  the  equillbrlim 

of  the  forces  exerted  on  the  rib  by  the  shell  is  to  be  ensured* 

Accordingly ,  solutions  obtained  in  this  section  are  for  u  ^  1* 

In  order  to  determine  U  .  V  ,  W  .  the  inverse  Fourier  transform. 

m  m  m 

defined  as 

f(n)  *  (2*)~^£  [f(a)/exp  inn]  da  ,  (18) 


is  applied  to  Eq.(16),  the  resulting  integrals  being  evaluated  by 
means  of  contour  integration  in  the  lover  half  of  the  complex  plane » 
faking  account  of  the  fact  that  is  an  odd  function  of  n  while 
and  are  even  functions.  Since  A  is  a  fourth  order  polynomial 
in  a2  with  real  coefficients,  it  can  be  written  as 


a  r  A  ,  (1)2  .  (3)2.  2  .  (1)2  (3)2,  f  4 

A  *  a.  [o  -  (a'  +  u '  )a  +  a A  &  }  (a  - 

o  mm  mm 


,  ,(4)2)a2 
m  m 


+  «(2>V4>2j 

m  a 


(s)2 

where  *  a  *  l,  2,  3,  4,  are  the  roots  of  the  polynomial,  such 
that 


(19) 


csmoua.  tkchhouxs  corporaiiok 


For  a  >  2,  it  follow  that 

2 

{V  V  V  ■  2P,RD"1^1  ta  ({«gn(n)  qj8>, 

■ 


where  ga  denotes  the  real  part,  and 

(Qy*>,  q{8),  li*(a){-flf  f2,  f3)  i(o  -  48>)M 

mam 


are  the  appropriate  residues  at  a  staple  pole,  the  numerically 
evaluated  roots  of  A,  for  the  corresponding  range  of  m  employed  In 
the  calculations,  having  been  found  to  be  distinct  and  different 
from  the  roots  of  f^,  fg,  and  f^. 


For  a  «  0.  f~/A  vanishes,  therefore,  *  0.  At  the  origin  of 

4  u  (2) 
the  conqplex  plane,  where  A  now  has  a  double  root  denoted  by  , 

fj/A  has  a  simple  pole,  while  f^/A  is  analytic,  so  that 

q52)-  W(v2  -*  -  1),  Q<2>  «>  0.  (23) 

0  0 


Furthermore,  in  the  lower  half  of  the  complex  plane,  A  has  the 
remaining  two  distinct  roots  and  a^,  different  from  those 
of  f^  and  f^,  and  the  corresponding  residues  can  be  evaluated  from 
Eq. (22)  with  v  *  0.  Then,U^  and  Wq  are  detezminad  from  Eq, (21)  with 

*  ■  0. 

Finally,  considering  a  given  number  N  of  pairs  of  arbitrary 
radial  line  loads  located  symmetrically  across  the  transverse  axis 
of  the  shell  at  n  *  J  t^/R,  J  »  1,  2,  ...  ,  N,  where  represents 
half  the  distance  between  the  j~th  pair  of  ribs,  the  solutiv  *or 
the  displacements  is 
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°b  “  Jo  X  "i  «*  'j 

*1  ■ 

+  sgn  n!J  exp  £j)]  coe  «•, 

vb  -  -2KflA  Jo  «j 

■«  * 


+  exp  5j)l  Psj«ln  «9» 


N  2 

«b  -  -2ED'1J£]  Jo  A  5' 

w»l  * 


♦  exp  ?'')]  PEjCoa  mB, 


r  (24) 


nj  *  »  n!j  *  n  ♦  t^/R  % 


(a) 


where  p  are  Interaction  load  amplitudes  for  thej-th  pair  of  radial 
line  loads.  This  solution  allows  non- uni fora  spacing  of  adjacent 
riba. 

Substituting  Eq,(24)  into  the  following  resultant  atresa- 
dlsplacement  relations13 


■bn"1*"1  V+"("b  +  Tby-V  ’ 


„-l 


NbB  *  DR  **b  +  Vb,8  +  vDb,n  +  ‘<Wb  +  Wb,08^' 


3,(1-  vJba'X.o  +  V^^V-W1* 


.-1 


’W  *  «l  -  *>DX  +  V  +  *(ub,e  +  Wlj 
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«M“KD<«b  +  «b,eb+vWb,nn)’  <«> 

^  -  (1  -  v)Kd(vm  -  Vn9). 

^-Ha-vkDCu.^-v^  +  aw,^), 

.  I 

the  corresponding  solution  for  the  stress  end  couple  resultants  is 

obtained*  These  results  ere  then  represented  in  the  p ,+  coordinate 

6  •  ■> 
system  by  transformation,  after  which,  expressions  for  and  N^ 

are  found  from  Eqs.(2)  and  (3)« 

4*  CYLINDRICAL  SHELL  UNDER  ARBITRARY  LOADING  ALONG  A  CIRCULAR  PENETRATION 

For  a  long,  unstlffened,  circular  cylindrical  shallow  thin  shell 

under  edge  loading  only,  the  stress  and  couple  resultants  can  be 
2 

expressed  in  terms  of  a  complex,  dimensionless  potential  function 
♦  “  , 

NC»  ■  bp"2Sb(**„  +  pV’  %  “  B^(**pp)’ 

%♦  ■  Nc,p  ■ 

Mcp  "  -*PotaI*»pp  +  p_2(*’^  + 

%  '  *^0— l'*’pp  +  0  +  P*’p^’ 

%  “  _i%p  “  (l  '  w)*Pop’2R«(p«.p<(  -  *>*>> 

B  spEt2/[l2(l  -  v2>f*  , 


where  Im  denotes  the  imaginary  part,  and 
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vV*  +  8i*V*,.vt  -  0,  MM  i  -Wc/P0, 

**<  )  2  (  ♦  <  ),w.  p2  3  *13(1  -  V2)!*/**, 


<a*> 


with  '  as  the  material-curvature  parameter*  Various  conditions 
concerning  the  applicability  of  the  shallow  shell  theory  to  the 
present  problem  have  been  recorded  elsewhere.*’ 

The  solution  of  Sq»(27)  in  the  p,4  coordinate  system, 
satisfying  the  conditions  of  biaxial  symmetry  and  of  decaying 
to  taro  m  p  4  •,  ia^ 


♦  “  d"Mo  Jo(a«  +  Ui)iMflc®'  2a*' 

\»l  *  (1  -  4a0>  WH2l  • 


(28) 


mhcre  u,i  are  Integers,  6  ~  is  a  Kronecker  delta f  A  ,  A*  are 

nu  m  *  ■ 

arbitrary  integration  const ante,  and  and  are  the  Bessel 
end  iankel  functions  of  the  first  kind  and  Integer  order  k, 

v 

ref  actively,  of  the  corplex,  dimensionless  argument  8p (2i); 

three  functions  are  evaluated  by  means  of  a  aeries  representation. 

Egressions  for  <)  and  N .  are  found  from  Eqe*  (2)  and  (3)  by 
cp  cpo 

si  ploying  Eqs.(26)  and  (2  8). 

In  order  to  determine  U  and  U  with  If  given  by  Eq.(2?), 

cp  e9  c  3 

tie  following  resultant  strees-dls placement  relations’2 


ucp5p  -  -  WcR_1«i«2 

+  %„>>  “  -  'V  "  HcR"1  COsZ  ♦* 

<=‘3<V»  -  V  +  °c*,P  “  2<*  +  -  WcR_1 


>  (29) 


sin  29 , 


at'j  utilized ,  the  solutions  being  subject  to  conditions  of  symmetry 
spj rifled  above. 

6 

Finally,  W  in  the  t^,0  coordinate  sycfeta  is 
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wc  “  ®lJo  X(cla;  -  c«v  **  -s«  <*» 

vfaerc 

C^Cn)  4*  iC^(n>  *  cos  t0i|R(2i)^  -  %iij  exp  %tnwf 

°«(B> 5  { Sir1}  J  Hn1}  «•  »* c~  a9  de  1£{«-oJ-  <31) 

noting  the  transformation  relations  in  Sq.(l);  the  integral 
applying  in  Be. (31)  le  evaluated  numerically. 


5*  BOUNDARY  CONDITIONS 


Considering  first  the  boundary  conditions  at  the  riba, 
since  the  total  solution  is  symmetric  in  n,  only  the  conditions 
at  n  -  i^/R,  i  -  1,  2,  ...  ,  N,  namely. 


W  -  w1  , 

13 

need  to  be  employed,  where 

w'  -  cu  cos  $  +  I'^E'D^J^ICl  -  I/AR'2)^ 

«K‘l 


(32) 


-  »“\d 


cos  n6 


(33) 


ia  the  otrward  radial  displacement  of  the  i-th  pair  of  ribs,  with 

C11  denoting  arbitrary  integration  constants  representing  rigid 

body  displacements.  From  Eqs. (9) ,<24) , (30) ,  (32),  and  (33)r_it 
7 

follows  that,  for  aP  1, 


N 


pnj  si hw' ' 2H1 ' 


»VW 


m0 


+  Jo(g™ia;  -  wv'Sc 


(34) 


vhere  B~i.  is  the.  Inverse  of 


5  <«ap  ^  +  «p  q,>i 

+  (DR’4<E’IR)[(x  -  I/AR,2)6  -  -  (m2  -  I)~2]6 


(35) 
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Csoi'  cLlf  Cji*  “*  5ji  *****  cw*  Cm*  5j>  *nd  5j* 
respectively,  evaluated  itj|-  i^/R,  and,  for  a  -  1, 

cu  -  m_lnlo<Ci»lA;  -  ClniV*  <*> 


Considering  next  the  reinforced  penetration,  the  loads 
exerted  on  the  ring  by  the  hydrostatic  pressure  and,  along  its 
external  and  internal  boundaries,  by  the  shell  and  t':s  membrane 

o 

covering  the  opening,  are  equipollent  to  the  following  system 
of  radial,  tangential, and  transverse  forces,  F^,  F^,  F^,  respec¬ 
tively,  and  torsional  sonant  M  acting  elong  p  •  p. 


P2p;S  «  *lPo'Pp  -  +  pWp  +  PoR'1(Np  ai*2  * 


The  equilibrium  of  these  forces  and  moment,  the  positive  directions 
of  which  can  be  inferred  from  Eq.(3  7)  and  Fig, 3,  is  essentially 
ensured,  since  it  can  be  shewn  that  Che  solution  in  Eq.(?8)  implies 
a  self-equilibrating  system  of  loads  along  the  boundary  of  the 
penetration. 

Denoting  by  u,  vp  and  w  the  displacements,  in  the  outward 
radial,  tangential,  and  upward  transverse  directions,  respectively. 


of  points  of  the  ring  on  p  ■  p,  the  boundary  conditions  et  the 
0 

penetration  are 


f\  “  +  **2(u  + 

-  %♦>  + 

V 

v\  *  ».pW)  -  EijCw, ^  +  s“^.****)  > 


2" "»S** 
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where,  epproxlaetely*, 

«  -  ♦  h(*)V,a  ♦  pjChl  .la2  +  , 

V  *  +  h^>w,4  +  »in  2*  • 

w  «  V  -  eV  . 

#P 


y  09) 


Shall  quantities  appearing  in  Eqs,(37)  and  (39)  are  evaluated 

at  p  »  p  ,  end,  in  the  course  of  the  solution,  h(+)  and  w  ^  are 

o  ,p 

replaced  approximately  by 

«»(♦)  -  h  -  JspV1  .in2  ♦,  h  s  h  (0),  w  •  *  . 

0  ,P  ,p 


6*  NUMERICAL  AND  EXPERIMENTAL  RESULTS 

The  above  analysis  has  bean  coded  and  numerical  results  for 

certain  configurations  within  the  range  of  applicability  of  the 

theory  and  for  which  experimental  results  are  presently  available 

have  been  generated  by  the  use  of  a  digital  computer. 

In  the  special  case  of  a  cylindrical  shell  stiffened  by 

multiple  ribs,  with  an  unreinforced  penetration,  a  comparison  of 

the  previous  analyses**’2  and  experimental  data  has  already  been 

reported. 22  These  analyses  have  also  been  confirmed  by  further 
14 

numerical  results  obtained  on  the  basis  of  an  energy  method  utilising 
a  finite  difference  representation.  The  following  account,  then, 
consists  of  similar  comparisons  between  the  results  of  the  present 
analysis  and  those  of  further  experiments  for  the  regaining  cases  of 
an  unstiffened  shell  with  a  reinforced  penetration  *  and  a  rib-stiffened 
shell  with  a  reinforced  penetration^,  respectively,  for  which  no  such 
evaluations  yet  exist.  The  corresponding  zone  of  influence  of  the 
penetration  and  the  state  of  stress  concentration  about  it  are  thereby 
delineated. 

Numerical  results  obtained  during  the  present  investigation 

for  the  {.articular  case  of  an  ur stiffened  cylindrical  shell  with 

a  reinforced  circular  penetration,  and  describing  the  variation 

of  the  dimensionless  hoop  and  longitudinal  stresses  with  p/p  at 

o 

the  outer  and  inner  surfaces  of  the  shell  for  $  »  0  and  w/2,  are 
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indicated  by  th*  solid  curves  In  Fig*  ,4-7  dwrt  tbs  experimental 
12 

results  for  the  same  cum  are  shown  by  the  dashed  curves.  The 

12 

latter  are  based  on  photoelasticity  tests  performed  on  an  epoxy 
(Hysol  4290)  shell  nodal,  with  a  ring  of  rectangular  cross  section 
sad  of  the  sans  material  around  the  opening  covered  with  a  plug, 
loaded  by  Internal  hydrostatic  pressure.  The  relevant  dimensions 
for  the  test  nodal  are  R  •  5-7/8  inches,  t  *  1/4  inch,  pQ  •  7/8  inch, 
h^  ■  1/2  inch,  and  e  ■  a  *  1/8  inch;  the  ring  has  a  fillet  with  a 
radius  of  1/16  Inch,  the  point  at  which  it  neats  the  shell  being 
indicated  by  a  vertical  dash  on  the  experimental  curves.  The  value 
of  E  is  taken  as  sero. 

Corresponding  results  for  the  general  case  of  a  rib-stiffened 
cylindrical  shell  with  a  reinforced  circular  penetration  are  given 
in  Figs .8-11  where  the  solid  curves  are  obtained  from  the  present 
analysis  and  the  dashed  curves  represent  the  experimental  data. 

The  latter  are  based  on  tests,  similar  to  those  mentioned  above, 
conducted10  on  a  model  of  the  same  material  and  dimensions,  except 
for  the  presence  of  multiple  ribs  of  rectangular  cross  section 
and  of  the  same  material,  with  adjacent  ribs  spaced  uniformly  a 
distance  21^  apart,  as  shown  in  Fig.l,  the  additional  relevant 
dimensions  being  l ^  •  1-1/2  inches,  t* »  5/16  inch,  and  e*  *  1-1/8 
inches,  where  t*  and  e'  are  the  width  and  the  depth  of  a  rib, 
respectively.  The  calculations  were  carried  out  for  the  case  N  »  4; 
the  centerlines  of  the  first  two  riba  only  are  shown  in  Figs. 8  and 
9,  with  the  width  of  the  ribs  indicated  by  a  pair  of  vertical  dashes 
on  the  experimental  curves. 

Finally,  the  calculated  values  of  the  principal  stress  concen¬ 
tration  factor  at  the  outer,  inner,  and  middle  surfaces  of  the 
shell,  for  the  two  cases  considered  here, are  shown  in  Fig. 12  for 
0  «£  $  ff/2  and  p/pQ  «  1,  with  the  curves  a  and  b  referring  respectively 
to  the  unstiffened  and  the  rib-etiffened  cylindrical  shell  with  a 
reinforced  circular  penetration. 
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7.  CONCLUSION 

From  the  preceding  graphical  comparisons  of  analytical  and 
sspsrlmsntsl  raault a,  it  can  be  concluded  that,  in  general,  a 
aatie factory  agreement  exist*.  Coaplete  agreement  is  not  expected, 
since  there  are  certain  features  of  the  test  models  causing 
localised  perturbations,  such  as  the  fillets  around  the  rin^  and 
the  ribs,  the  finite  width  of  the  ribs,  and  the  particular  plug 
covering  the  penetration,  that  are  not  accounted  for  in  the  theory. 
It  should  also  be  noted  that  sons  of  the  experimental  data  reported 
for  the  fillet  region  at  the  ring  are  based  on  extrapolations. 

It  may,  however,  be  desirable  to  refine  the  present  analysis 
by  considering  the  influence  of  the  circumferential  interaction 
between  the  shell  and  the  ribs. 
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